Are two kinds of dark matter seen in Galactic gamma rays? 
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Excesses of Galactic gamma rays in the 1-100 GeV region observed by EGRET and by WMAP 
have been interpreted as signals of dark matter (LSP) annihilation. We argue that the excess of TeV 
gamma rays from the Galactic center observed by H.E.S.S. signals a heavier dark matter component 
with a very small relic density. 

PACS numbers: 95.35.+d, 14.80.Ly, 98.35.Jk, 98.70.Vc 



I. INTRODUCTION 

Dark matter has so far manifested itself only through 
its gravitational effects - in fact on several scales in differ- 
ent ways. If dark matter is composed of stable, massive 
neutral particles having no or very weak interactions with 
baryonic and leptonic matter, WIMPs, it may be very dif- 
ficult to detect them unless they annihilate into photons, 
either directly or via quark-antiquark pairs or lepton- 
antilepton pairs. The popular Supersymmetric (SUSY) 
models Q offer several candidate particles, in particular 
the lightest SUSY particle (LSP) which must be stable. 

In our Galaxy dark matter is abundant, forming of 
a halo with some radial density distribution. Until 
particle physics experiments clarify the nature of dark 
matter, the best source of information may be annihi- 
lation gamma rays from the Galaxy center (GC), from 
the disk, the spiral arms, and from the halo. Indeed, 
several recent observations of gamma rays in different 
energy ranges |^ |^ |^ have been analyzed and inter- 
preted as signatures of LSP dark matter self-annihilation 



II. EGRET 

The diffuse Galactic 7-ray spectrum observed by the 
EGRET satellite extends from about 30 MeV to almost 
100 GeV [1 [13. Below 1 GeV the spectrum is fairly 
well described by conventional sources: energetic bary- 
onic interactions in the ISM resulting in secondary tt" 
decays, inverse Compton (IC) scattering of electrons on 
starlight and on CMB photons, and Bremsstrahlung from 
electrons |3 i, HH. Above 1 GeV the EGRET spec- 
trum shows a significant excess over the conventional 
sources, with the same spectral shape in 6 sky directions, 
as is expected for dark matter annihilation 8] through- 
out the Galaxy. Various explanations modifying one con- 
ventional component but not invoking dark matter have 
been tried (referenced and discussed in (ILj). but either 
the baryonic component has to originate from a harder 
proton spectrum, in conflict with local limits on antipro- 
ton and positron abundances, or the Compton scattering 
electrons have to be harder than what is locally observed, 
or one has to adjust the absolute normalizations of all the 



three components ad hoc. 

Invoking the presence of dark matter, a hard electron 
spectrum could be obtained from LSP annihilations xXi 
mainly into hh quark pairs. Synchrotron radiation from 
such ultrarelativistic electrons, as evaluated with the help 
of the DarkSusy program |l2l | and the Pythia program 
[T^j l. turns out to reproduce the shape of the EGRET 
excess spectrum well for LSP masses in the range = 
50-90 GeV jlj. 

However, to obtain the absolute normalization in each 
of the 6 directions one would need a model for the triax- 
ial halo profile of dark matter densities. Although there 
exist phenomenologi cal p rofiles derived by averaging over 
hundreds of galaxies the average profile is not 

necessarily a good description of our Galaxy. Inversely, 
one can use the observed excess spectrum to trace the 
profile: a spherical WIMP halo radiates isotropically, the 
higher star density in the disk testifies to an enhanced 
WIMP density there, and the GC density may be more 
or less cuspy. De Boer et al. @ propose a disk model 
boosted with two radial rings of enhanced density, the 
inner one fitted at radius 4.3 kpc and the outer one at 
14 ± 2.1 kpc. It is intriguing to note that this model is 
supported by the observation of visible matter forming a 
ring of stars at 14 — 18 kpc with enhanced density. This 
two ring model also fits the rotation curves well, includ- 
ing the mysterious change of sign of the slope at 9.4 kpc 
(for references see jsj). 

Thus the EGRET excess appears to be well explained 
by 7-rays from the self-annihilation of LSPs with a 
relic density in consistency with the first year result of 
WMAP [itI, = 0-23 ± 0.04, and with a total annihi- 
lation cross-section of (ctaw) ^ 2 ■ 10~^^ cm^/s . 

III. WMAP 

Let us now turn to another window. In a care- 
ful analysis of Galactic foreground components in the 
WMAP data, Finkbeiner has noted an excess emis- 
sion from the inner Galaxy, distributed with radial sym- 
metry within 20° — 30° around the GC. He argues that 
this may represent inverse Compton scattered starlight 
on diffusion hardened electrons, produced by annihila- 
tion of ~ 100 GeV WIMPs into e+e" pairs. 

To analyze this radiation Finkbeiner R| assumed a 
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smooth spherical density distribution with a truncated 
NFW profile 0| in the inner 1-2 kpc Galaxy, and 
evolved the electron spectrum through diffusion and IC 
energy loss. 

The resulting 7-ray spectrum appears to have a shape 
rather similar to that of the EGRET excess (an exact 
comparison is rendered difficult by the different units 
used in 7] and 8]), thus supporting the interpretation 
of LSP annihilation. The two analyses differ clearly in 
their density profiles in the galactic plane, and also in the 
annihilation mechanism: the power is assumed either to 
go primarily into bb Q or entirely into ultrarelativistic 
e+e~ pairs One can still conclude that the two stud- 
ies observe the annihilation of the same particle x- 

IV. H.E.S.S. 

Aharonian et al. 6] have announced the detection of 
very high energy (VHE) 7-rays from a point-like source 
coincident within 1' of Sagittarius A* in the GC, obtained 
with the H.E.S.S. array of Cherenkov telescopes. The 
observed energy spectrum ranges from 0.3 TeV to 10 TeV 
and is rather flat. 

Although various conventional sources of VHE 7-rays 
may still be plausible explanations we shall take 
the attitude that the source is annihilation of a heavier 
WIMP, ^, which cannot be the source of the much softer 
X annihilation 7-rays seen by EGRET and WMAP. If x is 
the LSP, as has been argued above and in J,,^], the much 
heavier ^ would have frozen out of thermal equilibrium 
earlier during the radiation-dominated epoch. 

Let us take the ^ mass to be toj, the relic density 
parameter f2j and the average annihilation cross-section 
(aAv). The ^ fall out of thermal equilibrium at temper- 
ature Teq <^ mj, and they annihilate until their number 
density has been reduced to give {aAv)ni^ « H. The 
properties of the ^ are constrained by the relations 

oc {aAv)^^ {<^Av) oc m'^^ . (1) 

Using the WMAP value for the x relic density fi^ , a value 
for can be derived. Now, however, there is not room 
for a high ^ density as well as the high x density, so we 
must have ^ Cl^. In fact, the order of magnitude 



is set by the error on the ft^ measurement '17||, thus 

< 0.04. 

This relic density has already been estimated by Silk 
for a typical NFW profile and an inferred cross- 
section of 10 pb at the upper end of the H.E.S.S. energy 
spectrum, 10 TeV, with the result fi^ w 0.03. This is 
exactly the value needed, although it has met with great 
difficulties in previous analyses 6, 9, lOj which treated 
the H.E.S.S. detection as LSP annihilation, not recogniz- 
ing the role of the EGRET and WMAP excesses. 

Two- WIMP models that have lately been discussed 
are, for instance, by Boehm and Fayet who con- 
sider theories displaying an underlying N = 2 extended 
supersymmetry, and by Hooper and March- Russell |[l8l |. 
who take a stable dark matter candidate from the mes- 
senger sector of gauge mediated supersymmetry breaking 
models. With a mass of 20 to 30 TeV they can indeed 
generate the spectrum seen by HESS with the WMAP 
relic density. They probably possess enough freedom also 
to match a very small relic density, as we have advocated 
here. 



V. DISCUSSION 

This two- WIMP conjecture has several interesting con- 
sequences. Knowing the approximate masses and anni- 
hilation cross-sections of two neutral SUSY particles is a 
considerable new input to SUSY spectroscopy. 

Having two dark matter particles with so different 
masses changes the picture of galaxy density profiles. 
The X need not have much of a central cusp, whereas 
the ^ probably dominates the innermost region. This 
would feedback to the EGRET analysis which needed an 
inner ring of dark matter Q, as well as to the WMAP 
analysis which used a truncated NFW profile . 

The picture is considerably broadened if one also con- 
siders schemes of coannihilation. Also, some of the 
charged ^ annihilation products may perhaps diffuse out 
of the GC and spill lower energy photons into the upper 
ends of the EGRET and WMAP spectra, thus affecting 
the analyses at hand. There is clearly need for a simul- 
taneous study of the EGRET excess and the WMAP ex- 
cess. 
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